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Measurements of on-shell ZZ production are described, using data from the 
DELPHI experiment at LEP in e^e~ collisions at centre-of-mass energies be- 
tween 183 and 209 GeV, corresponding to an integrated luminosity of about 
665 pb^^. Results obtained in each of the final states gggg, vvqq, fi'^fi~qq, 
e'^e~qq, r^r^qq, and vvl'^l' (with / = e,/i) are presented. The 

measured production cross-sections are consistent with the Standard Model ex- 
pectations. These results update and supersede those already published at 183 
and 189 GeV. 
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1 Introduction 

The first evidence [1] for doubly resonant production of Z bosons was observed during 
1997, when the LEP-2 accelerator reached a centre-of-mass energy near 183 GeV, corre- 
sponding to the threshold for this channel. In this article we present measurements of the 
production cross-section, using data collected by the DELPHI experiment in 1997-2000, 
at centre-of-mass energies up to 209 GeV, corresponding to an integrated luminosity of 
about 665 pb~^. 

There are several motivations for studying this channel. Firstly, it enables a check of 
the Standard Model (SM) prediction. Both the cross-section and angular distribution of 
the produced Z bosons are sensitive to contributions from new physics beyond the SM. 
Several mechanisms for anomalous production [2,3] can hence be constrained directly. 
For instance, limits on neutral triple gauge boson couplings [4], which are forbidden in 
the SM, can be set [5]. Secondly, the ZZ production process is an almost irreducible 
background to the Higgs search at LEP, particularly when the mass of the Higgs boson 
is close to that of the Z [6]. It provides an environment similar to that corresponding 
to a possible Higgs signal, in the main final state topologies which arc analysed, both in 
terms of experimental signature and rate. In this context, it allows to test the techniques 
used in the Higgs analyses on an existing process. 

In what follows, the data sets and simulations used are described and the signal def- 
inition adopted is discussed. The event selections are presented for each of the seven 
sub-channels which were analysed: qqqq^ vvqq^ ^~^^~qq, e~^e~qq, T~^T~qq, and 
vvl^l^ (with / = e, fj,). Results are given in the form of a comparison of the numbers of ob- 
served and predicted selected events, together with an evaluation of the main systematic 
effects. Finally, combinations of sub-channel results into overall ZZ cross-sections at each 
centre-of-mass energy are described, as well as combinations performed over all centre- 
of-mass energies, both separately in each sub-channel, and for all channels together. The 
resulting measurements are compared with SM expectations. 

These results update and supersede those already published at 183 and 189 GeV [7], 
by including the T'^T~qq sub-channel and using improved methods in several of the other 
channels. Measurements of on-shell ZZ production published by the three other LEP 
collaborations can be found in [1,8-10]. 

2 Data samples and simulation 

The data samples used were collected by DELPHI in the years 1997-2000. The cor- 
responding centre-of-mass energies and integrated luminosities which were analysed are 
given in Table 1. A detailed description of the detector and a review of its performance 
can be found in [11,12]. For LEP-2 operation, the vertex detector was upgraded [13] and 
a set of scintillator counters was added to veto photons in blind regions of the electro- 
magnetic calorimetry, at polar angles 9 ~ 40°, 90° and 140°. 

Simulated events were produced with the DELPHI simulation program DELSIM [12] and 
were then passed through the same reconstruction and analysis chain as the data. The 
generation of processes leading to four-fermion final states was done with EXCALIBUR [14] , 
relying on JETSET 7.4 [15] for quark fragmentation. GRC4F [16] was used as a comple- 
mentary generator for four-fermion final states resulting from single-resonant Weve and 
(Z/7*)e+e^ processes when the spectator electron was close to the beam direction. Two- 
fermion processes e"'"e~ — >• 55(7) were generated using PYTHIA [15], e^e'^ ji^ ji^ij) and 
e+e- T+T-{'y) with KORALZ [17], and e+e' e+e-{-f) with BHWIDE [18]. Two-photon 
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interactions were generated using TWOGAM [19] and BDK [20]. The sizes of the simulated 
event samples were typically at least 100 times larger than expected in real data. 

One sector (corresponding to 1/12) of the main tracking device, the Time Projection 
Chamber, was inactive for the last quarter of the 2000 data sample. The corresponding 
small change of analysis sensitivity from this period was studied with specially made 
simulations and taken into account in the extraction of the cross-sections. 



Year 




Integrated 




[GeV] 


luminosity [pb~^] 


1997 


182.6 


54.0 


1998 


188.6 


158.1 


1999 


191.6 


25.8 


1999 


195.5 


76.9 


1999 


199.5 


84.3 


1999 


201.6 


41.1 


2000 


< 205.5 


83.3 


2000 


> 205.5 


141.8 


Total 




665.3 



Table 1: Centre-of-mass energies and integrated luminosities of the analysed data. During 
the year 2000, the energies reached were in the range 202-209 GeV, clustered mainly 
around 205 and 207 GeV. 



3 Signal definition in the simulation 




Figure 1: The Feynman graphs for on-shell ZZ production (referred to as the NC02 
graphs) 



The region of phase-space at high di-fermion masses which characterizes ZZ produc- 
tion must be isolated to measure the corresponding cross-section. In this region, be- 
sides the dominant tree-level doubly-resonant ZZ production graphs shown in Figure 1 
(referred to as the NC02 graphs), there arc also contributions from other four-fcrmion 
processes, which can lead to identical final states^. In order to interpret the measure- 
ments in terms of NC02 graphs, the ZZ signal was defined in each of the sub-channels 

^Such contributions arise from e+e^ (Z/7*)7* processes (when 7* virtualities arc close to the Z boson mass), from 
(Z/7*)e+e~ processes (in cases of final states with electrons), and from processes involving W bosons (in cases of uudd, 
cess, and vivil'^l~ final states). 
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by pre-selecting a sub-set of events within the four-fermion simulation, having both a 
high purity and efficiency in terms of the relative contribution from the dominant NC02 
graphs. In the analyses of l^l~qq (with I — e,ii,T) final states and in one of the two 
analyses of the uuqq final state (the probabilistic selection) described in Sections 5-7, 
the signal definition used was based on cuts on the generated boson masses, which were 
required to be within 10 GeV/c^ of the nominal Z mass [7]. In the other sub-channels, 
the signal was defined as the sub-set of events satisfying: 

^^7^ > 0.5, 
\Mau\ 

where A1nco2 and Maii are the values of the matrix elements, computed using the four- 
momenta of the generated fermions, for NC02 and for all four-fermion graphs, respectively. 
Corrections were necessary to relate the rates of such events to pure NC02 cross-sections^, 
in order to take into account the efficiency of the signal-defining selection cut, the residual 
contamination from the non-NC02 processes, and interference effects. However, because 
NC02 and non-NC02 contributions are naturally well separated in phase-space, the mag- 
nitudes of these corrections were typically less than a few percent. Simulated four-fermion 
events which did not satisfy the signal- defining cut were considered as background in the 
comparison with data. 



4 Four jet channel 

The qqqq decay mode represents 48.9% of the expected ZZ final states. It results 
typically in four or more well separated jets of particles. The dominant backgrounds 
originate from fully hadronic WW final states and from qq{'y) processes with hard final 
state gluon radiation, both of which can lead to similar multi-jet topologies. The main 
ingredients required to isolate the signal were di-jet mass reconstruction, topological 
information quantifying the jet separation and kinematics and (for bbqq decay modes) 
the 6-tagging of the jets. A probabilistic method based on likelihood ratio products was 
developed to combine in an optimal way the information and compute an event-by-event 
measure of the compatibility with the ZZ hypothesis, taking into account the flavour 
content (with or without 6-quarks) of the signal. 

4.1 Event pre-selection 

A pre-selection was applied to select fully hadronic events. Firstly, events were required 
to contain at least 18 reconstructed tracks of charged particles and to have more than 
69% of the available centre-of-mass energy observed in the detector. To reduce the 
contamination from 55(7) processes with energetic photons emitted in the beam pipe, 
the effective centre-of-mass energy of the event, y/s', computed as described in [22], was 
required to exceed 80% of ^/s. Particles were then clustered into jets using the DURHAM 
clustering algorithm [23] (with ycut=0.001) and events were selected if at least four jets 
were reconstructed. To reduce backgrounds from processes with energetic and isolated 
leptons or photons, jets were required to contain at least four particles and have an 
invariant mass above 1 GeV/c^. The efficiency of these cuts was about 88% and 80% for 
the fully hadronic ZZ signal and WW background, respectively, while about 3% of events 

^Destructive interference from Fermi correlations occurs in the case of final states with four identical fermions [21]. 
Although modelled by the four-fermion simulation, the resulting reductions in overall cross-section, smaller than 0.5%, 
were neglected in the calculation of these NC02 cross-sections. 
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from qq{'j) processes remained. At this stage, backgrounds from all other processes were 
negligible. 

4.2 Probabilistic selection 

The method developed enabled the computation of the relative probability for any 
pre-selected event to originate from each of three SM processes, WW, ZZ or 55(7). For 
this purpose a calculation based on likelihood ratio products was used, based on the 
main features characterising the ZZ signal: reconstructed di-jet masses, jet kinematics 
and separation and possible presence of 6-hadrons in the jets. These three ingredients 
and their combination are described below. The computation combined measurements of 
the events with the theoretical expectations corresponding to each hypothesis, described 
analytically where possible in order to reduce the dcpcndance on detailed simulations. 
Throughout the analysis events with four jets and events with five or more jets^ were 
treated separately. 

Reconstructed mass distributions can be used to characterise the different hypotheses. 
The expected shapes are well defined, but in reconstructing events ambiguities arise, since 
in a four or five-jet event, there are three or ten possible jet-jet pairings, respectively. In 
addition, the quality of the energy flow reconstruction can vary significantly from event 
to event. Following the scheme described in [24], these considerations were taken into 
account by simultaneously evaluating the compatibility of each event with each hypothesis 
considered, WW, ZZ and qq{'~f), from the mass information. For this purpose, probability 
distributions for the different mass combinations were first expressed in analytical form, 
as products of Breit-Wigner distributions and phase-space factors to describe the case of 
correctly paired mass combinations in genuine WW and ZZ events, and as flat spectra 
for all the other cases. Then, for each reconstructed event, the four-momenta of the 
jets and their estimated errors were used in kinematic fits [25], requiring four- momentum 
conservation and equality of the masses of the two jet systems'' with the pairs of values 
to be tested. The two-dimensional distributions of the probabilities obtained from 
these fits as a function of the pair of test values were then multiphed by the expected 
probability distributions and the results integrated to quantify the compatibility of each 
pairing with each of the WW, ZZ and qq{j) hypotheses. To compute global event-level 
compatibilities with each hypothesis from the mass information, these results were then 
summed over the possible pairings (assuming equal a priori weights) and normalised. 
As an illustration, the result obtained for the WW hypothesis is shown in the left-hand 
plot of Figure 2, for the pre-selected events with four jets. In this method, the inherent 
ambiguity in the pairing of the jets, one of the main difficulties in experimental analyses 
of multi-jet events, did not need to be resolved. 

Specific features of the jet kinematics and separation in events arising from four- 
fermion and ^'^(7) processes with hard final state gluon radiation leading to a four-jet 
event were also used. Because gluon emission off quarks has infrared and collinear di- 
vergences, these types of events tend to have cigar-like shapes, while the four-fermion 
processes are more spherical. This difference was exploited to construct a topological 
variable: Dp^r = El9l^/E^, where Ei {E2) is the leading (subleading) jet energy and 
9i {62) is the smallest (next to smallest) opening angle between two jets. The correspond- 
ing ratio of probability density functions for these two event types was parametrized as 
a function of this variable using simulated data for four and five jet events separately. 

''Events with more than five jets were forced into a five-jet configuration. 

■'Two di-jets in the case of a four-jet event, and both a di-jet and a three-jet system in the case of a five-jet event. 
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Figure 2: The left-hand plot shows the distribution of the variable quantifying for each 
event the compatibility with the WW hypothesis, computed using only mass information. 
The right-hand plot shows the distribution of the event-topology variable Dpur (for events 
with Dpur/lO greater than 7.5 GeV^^'^rad'^/^). Both distributions are for pre-selected 
events with four jets (corresponding to about 80% of the total sample). The dots are the 
data taken in 1997-2000 and the histograms represent the simulation. 



The distribution of the Dpur variable for four-jet events is shown in the right-hand plot 
of Figure 2. 

The presence of 6-hadrons in the jets was identified using a dedicated algorithm [26] 
which exploited their characteristically long lifetime and high mass. This was made 
possible by the capabilities of the vertex detector [13], in which tracks originating from 
the displaced 6-hadron vertices could be resolved. The combined 6-tagging variable which 
was constructed is shown in Figure 3 for jets in hadronic Z decays (left), and in the pre- 
selected fully hadronic high-energy event sample (right). It was particularly powerful in 
reducing the contamination from WW events, since the corresponding final states hardly 
ever contain a 6-quark, while 38.7% of fully hadronic ZZ events have at least two 6-quarks. 
In order to use the information optimally, the ratio of the distributions of this variable 
for jets originating from the fragmentation of b and non-6 quarks was parametrized in 
three different angular regions, using simulated Z decays, to take into account the polar 
angle dependence of the detector resolution. Using these parametrizations and the values 
measured for each jet, relative b and non-6 probabilities could be assigned to each jet. 

Finally a variable quantifying the compatibility of each event with the ZZ hypothe- 
sis was constructed using a method based on likelihood ratio products to combine the 
information from the different ingredients described. Its distribution is shown in Figure 
4, including the data from all centre-of-mass energies. The different energies could be 
added without diluting the information because the distributions of the purity of ZZ 
events with respect to the constructed variable were the same. 
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Figure 3: The left-hand plot shows the distribution of the combined 6-tagging variable for 
jets in Z decays collected in special calibration runs during the high-energy data taking. 
The right-hand plot shows the distribution of the combined 6-tagging probability per jet 
for all reconstructed jets after applying the pre- select ion. The dots are the data taken in 
1997-2000 and the histograms represent the simulation. 



4.3 Results 

As an illustration, the observed and predicted numbers of events selected after a cut 
on the ZZ probability which maximises the product of efficiency and purity are shown 
in Table 2. There was an overall agreement between data and simulation within the 
statistical fluctuations. At each centre-of-mass energy, a measurement of the production 
cross-section was obtained from a binned maximum likelihood fit to the distributions of 
the ZZ probability, with the ZZ signal contribution as the only free parameter. The 
results obtained and the combinations performed between energies and with other chan- 
nels to derive global values for the NC02 cross-section are described in Section 10. The 
optimal region of the ZZ probability distribution to be used in the fit was determined by 
minimising the combined statistical and systematic uncertainty, evaluated using the full 
LEP-2 data sample. This resulted in a lower cut on the ZZ probability of 0.25. At this 
level, the number of events from the WW ^ QQil) and ZZ processes were comparable. 

4.4 Systematic uncertainties 

The main source of systematic error in the selection of ZZ — > qqqq events was from the 
limited precision available in the modelling of the multi-jet qq{'y) processes (particularly 
those involving 6-quarks) which comprised the main background in the most signal-like 
regions relevant to the extraction of ZZ cross-sections. Theoretical uncertainties in the 
predictions and biases from the generator treatments were studied at y/s = Mz, where 
Mz is the mass of the Z boson, and extrapolated to LEP-2 energies [27]. The most 
relevant aspects for the measurement of the qqqq channel are summarised below. 
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Figure 4: Distribution of the combined ZZ probability. The dots are the data taken in 
1997-2000 and the histograms represent the simulation. 



It was found that PYTHIA [15] underestimated the inclusive four-jet rate^ by typically 
about 10% at ^/s = Mz- The sign and magnitude of this discrepancy were also confirmed 
by comparisons at LEP-2 energies. After correcting for 1/2 of the discrepancy, ±1/2 of 
it was used as a conservative estimate of the uncertainty in the prediction of this rate. 

The probability of secondary c and 6-quark pair production through gluon splitting 
processes was found to be underestimated in PYTHIA by factors of about 1.5 and 2, re- 
spectively, compared to other calculations and dedicated measurements. An approximate 
correction of this deficit was achieved in the analysis by reweighting simulated events 
containing gluons splitting into heavy quarks with these factors. Since the gluon split- 
ting process remains poorly known both theoretically and experimentally, and since the 
reweighting procedure applied provided only a rough correction, a ±50% relative error 
on these probabilities was assumed. 

The magnitude of the reduction in gluon radiation off 6-quarks relative to other flavours 
(arising from their heavier mass) was shown to be overestimated in the version of PYTHIA 
used, compared to analytic calculations and to dedicated measurements, typically by as 
much as the theoretical uncertainty in these calculations. An additional ±4% relative 
error in the four-jet rate from gg(7) processes was conservatively assumed to cover both 
this uncertainty and bias. 

^The agreement between data and simulation for general hadronic event properties and shapes depends on the strategy 
applied in the tuning [28] of the generator parameters [27] . 
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ZZ 


mm 








Integrated 


Selection 


Predicted 


Predicted 


Selected 


[GeV] 


1 • • J rill 

luminosity [pb J 


efficiency 


background 


total MC 


data 


183 


54.7 


0.13 


0.92 


1.84 


2 


189 


158.0 


0.27 


17.88 


31.13 


29 


19z 


25.9 


O.o4 


5.31 


8.64 


11 


196 


76.9 


0.38 


21.16 


34.01 


46 


200 


84.3 


0.39 


26.64 


42.09 


36 


202 


41.1 


0.42 


14.24 


22.55 


26 


205 


82.0 


0.41 


27.45 


43.73 


45 


207 


142.2 


0.42 


52.58 


84.13 


78 


Total 


665.1 


0.38 


166.18 


268.11 


273 



Table 2: Integrated luminosities, selection efficiencies and number of observed and ex- 
pected events, after a cut on the probability of the ZZ hypothesis maximising the product 
of efficiency and purity. The measured ZZ production cross-sections were extracted at 
each centre-of-mass energy by fitting the distributions of this probability (see the text). 

The impact of propagating these uncertainties in the background level to the fitted 
ZZ — > qqqq cross-sections is shown in Table 3, using the full LEP-2 data sample. The 
effect from uncertainties in selection efficiencies related to the 6-tagging procedure is also 
shown, as well as the effect from varying the WW cross-section within an estimated 
uncertainty of ±2% [21]. The estimate of the error from 6-tagging was obtained by 
propagating the errors evaluated in [29] to the computed h and non-6 jet probabilities. All 
errors were assumed to be fully correlated between the centre-of-mass energies analysed. 





Inclusive 


Gluon splitting 


Gluon radiation 


6-tagging 


WW 


[GeV] 


four-jet rate 


into heavy quarks 


off' 6-quarks 


procedure 


cross-section 


183 - 208 


±2.3 % 


±2.3 % 


±1.8 % 


±1.5 % 


±0.9 % 



Table 3: Expected relative uncertainties on the fitted ZZ qqqq cross-section from 
systematic uncertainties in the predicted background and selection efficiencies. 



5 Jets and missing energy 

The vVqq decay mode represents 28.0% of the ZZ final states. The event topology 

is characterised by a pair of acoplanar jets (acollinear in the transverse plane to the 
beam) with visible and recoil masses compatible with the Z mass. Because two energetic 
neutrinos escape detection in this channel, efficient and reliable energy flow reconstruction 
is essential to select the signal. The most difficult physical backgrounds arise from single- 
resonant WeVe processes, from WW processes where one W decays into tPt, and from 
qq events, which may or may not be accompanied by energetic isolated photons escaping 
detection, in which one or both jets are badly reconstructed. 
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After a common pre-selection, two independent analyses were carried out. Each used a 
separate selection of discriminating variables to construct a combined estimator, employ- 
ing different, complementary schemes. The first analysis used a non-linear discriminant 
method (the so-called Iterative Discriminant Analysis (IDA) [30]) in which the cross- 
products between the variables were included, and where the correlations - important in 
this channel - can be treated optimally. The second analysis used a more conventional 
probabilistic method based on likelihood ratio products, in which the non-Gaussian pro- 
files of the variables are taken into account. Results from both analyses were obtained 
and compared. The IDA analysis had slightly better performance, and was chosen to 
derive the combined values of the NC02 cross-section described in Section 10. It was also 
used to study the propagation of systematic uncertainties. The probabilistic analysis 
served as a cross-check. 

5.1 Event pre-selection 

Selecting carefully the measured particles from which event variables are calculated 
is important to obtain efficient and reliable energy flow reconstruction. The definition 
of neutral particles was studied specially in this context, using both energy clusters in 
the calorimeters not associated to charged particle tracks and reconstructed vertices of 
photon conversions, of interactions of neutral hadrons and of decays of neutral particles 
in the tracking volume. Beam related backgrounds were also rejected by requiring at least 
two charged particles with impact parameter (with respect to the fitted primary vertex) 
less than 1 mm in the transverse plane and less than 3 mm along the beam axis, and 
with transverse momentum, P^, greater than 2 GeV/c. A loose hadronic pre-selection 
was then applied, requiring at least eight charged particles, total charged energy greater 
than 0.16\/i, transverse energy (defined as EP^, where the sum is over all particles) 
greater than 0.15^/s, and the sum of the component of momentum of all particles along 
the thrust axis greater than 0.25^^. Finally, events with an electromagnetic shower of 
energy exceeding 0.45^^ were rejected. This removed about 97% of the background from 
two-photon and Bhabha processes. 

The particles of the events were forced into two-jet configurations using the DURHAM 
jet algorithm. To reject events coming from qq{'-f) processes with energetic photons 
emitted in the beam pipe, the effective centre-of-mass energy of the event, y/s', computed 
from the jet directions, was required to be greater than 115 GeV when the polar angle 
of the total momentum of detected particles was less than 40° or larger than 140°. To 
reduce the contamination from ^'^(7) processes with energetic photons in the detector 
acceptance, events were rejected if their total electromagnetic energy within 30° of the 
beam axis was greater than 0.16^^, or if the total energy in the luminosity monitor 
was greater than O.OSy^. To reject qq{'y) processes with energetic photons emitted in 
blind regions of the electromagnetic calorimetry (at polar angles 6 ~ 40°, 90° and 140°) 
signals from the set of dedicated scintillator counters were used. To reduce the two- 
fermion background without significant initial state radiation as well as four-fermion 
backgrounds without missing energy, the effective centre-of-mass energy of the event had 
to satisfy y/s' < 0.96^/s. To reduce backgrounds from two-fermion events with jets 
pointing to the insensitive regions of the electromagnetic calorimeters, where the energy 
flow was less precise, events were rejected if the jets had polar angles within the ranges 
35 — 45° or 135 — 145°, unless the acoplanarity of the event was greater than 10°. This 
removed about 88% of the total qq{'y) background. 
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The contamination from WW processes where one W decayed leptonically was reduced 
by requiring that the energy of the most energetic particle of the event be less than Q.2^/s. 
To improve the rejection when the W decayed into ti/t, the events were required to have 
no charged particle with a transverse momentum with respect to its jet direction greater 
than 10 GeV/c. This removed about 66% of the WW background. 

At the end of the event pre-selection, the efficiency for the ZZ — > i/i/qq signal was 
about 77%. 

5.2 IDA selection 

After the common pre-selection described in Section 5.1, it was additionally re- 
quired that the momentum of the most isolated particle in the event be in the range 
0.01 < p < 0.20^/s, and that the visible mass be smaller than ^/s. A combined 
discriminant variable was then constructed using the Iterative Discriminant Analysis 
program (IDA) [30] to calculate a second order polynomial from twelve event variables, 
selected according to their discriminating power and independence ^: 

• The logarithm of the transverse momentum of the event with respect to the beam 
axis, for values greater than (see Figure 5 upper left), 

• The visible energy of the event normalised to ^/s, for values lower than 1 (see Figure 5 
upper right), 

• The transverse energy of the event normalised to ^/s, for values in the range 0.15 — 
0.60, 

• The minimum polar angle defining a cone in the positive and negative beam direc- 
tions containing 6% of the total visible energy, 

• The difference of the jet energies normalised to y/s, for values lower than 0.35, 

• The sum of energies of particles having angles with respect to the most isolated 
particle larger than 5° and smaller than 60° (25° if the momentum of the most 
isolated particle exceeded 5 GeV/c), normalised to ^/s, for values lower than 0.15, 

• The effective ccntrc-of-mass energy V^, 

• The logarithm of the product of the acoplanarity with the inter-jet angle, for values 
in the range 0.0 — 2.8 (see Figure 5 lower left), 

• The logarithm of the acoUinearity of the two jets, 

• The longitudinal momentum with respect to the thrust axis, normalised to y/s, for 
values lower than 1, 

• The missing mass normalised to ^/s, computed with the constraint that the visible 
mass equal Mz (see Figure 5 lower right), 

• The logarithm of the largest transverse momentum of any particle with respect to 
its jet direction. 

The observed distributions of these twelve variables were in good agreement with the 
predictions from simulation at each centre-of-mass energy analysed. As an illustration, 
the distributions of four of them are shown in Figure 5 after all pre-selection cuts, for 
^/s = 199.5 GeV. The IDA method for variable combination was applied in two succes- 
sive optimisations. Following the first one, an enriched sub-sample of events was selected, 
with a 95% relative signal efficiency. These events were then used in the second optimi- 
sation. The final combined discriminant variable obtained is shown in the left-hand plot 
of Figure 6, including the data from all centre-of-mass energies. 

''To concentrate on the signal region when optimizing the second-order discriminant function, and to avoid long non- 
Gaussian tails, additional very loose pre-selection cuts were applied to some of the variables, as indicated in the description. 
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Figure 5: Distributions of four variables used in the IDA analysis stream shown after 
all pre-selection cuts (see text): logarithm of the transverse momentum (in GeV/c) with 
respect to the beam axis (upper left), ratio of the visible energy to the centre-of-mass 
energy (upper right), logarithm of the product of the inter-jet angle (in degrees) and 
the acoplanarity (lower left), ratio of the missing mass (computed using the constraint 
that the visible mass be compatible with Mz) to the centre-of-mass energy (lower right). 
The dots are the data taken in 1999 at a centre-of-mass energy ^/s = 199.5 GeV. The 
histograms are the simulation prediction. 



5.3 Probabilistic selection 

After the common pre-selection described in Section 5.1, it was additionally required 
that the total visible energy be lower than 0.70a/s. A probabilistic method based on 
likelihood ratio products was then used to combine nine discriminating variables: 

• The energy of the most energetic jet, 

• The cosine of the acoplanarity, 

• The effective centre-of-mass energy y/s', normalised to ^/s, 

• The missing momentum, 

• The cosine of the polar angle of the missing momentum, 

• The largest transverse momentum of any particle with respect to its jet direction, 

• The minimum charged multiplicity of the jets^, 

^In this case events were not forced into two-jet configurations and the DURHAM jet clustering algorithm was used 
with ycut = 0.005. 
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Figure 6: Distribution of the combined discriminating variable for the vvqq analysis in the 
IDA (left) and probabilistic (right) selections. The dots are the data taken in 1997-2000 
and the histograms represent the simulation. The distributions of the purity of ZZ events 
as a function of these two variables were almost identical for all centre-of-mass energies, 
allowing to add the histograms of the different energies with minimal dilution of the 
information. 



• The visible mass, computed with the constraint that the missing mass equal M^, 

• The output of a veto algorithm based on the response of the dedicated scintillator 
counters installed at polar angles 9 ^ 40°, 90° and 140°. 

There was good agreement between the data and the predictions from the simulation 
for each of these variables at each centre-of-mass energy. The probability density func- 
tions used to construct the combined discriminant variable were determined from the 
simulation, separately at each energy. The distribution of this final variable is shown in 
the right-hand plot of Figure 6, including the data from all centre-of-mass energies. 

5.4 Results and comparison 

The results from the two selections are presented in Table 4, where the observed and 
predicted numbers of events, selected after cuts on each discriminant variable maximising 
the product of the efficiency and purity, are shown^. The two methods gave compatible 
results, but the IDA selection performed slightly better. There was overall agreement 
between data and simulation within statistical errors in both cases. A study of the 
overlap between the two selections, performed with the data collected in 2000, showed 
that the proportion of common events was about 60% in real data: 55% in the background 
and 85% in the signal, according to the simulation. The IDA selection was used as the 
main analysis, and the probabilistic as cross-check. 

At each centre-of-mass energy, a measurement of the production cross-section was 
obtained from a binned maximum likelihood fit to the distribution of the final IDA 

**The luminosities used to analyse the ZZ — > vuqq channel were slightly lower than for other channels, because additional 
criteria on the operation of the apparatus were used. 
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ZZ vvqq 
(IDA selection) 




Integrated 


Selection 


Predicted 


Predicted 


Selected 


[GeV] 


luminosity [pb^"*^] 


efficiency 


background 


total MC 


data 


183 


52.4 


0.31 


1.6 


3.1 


2 


189 


152.8 


0.49 


14.8 


28.3 


22 


192 


24.9 


0.56 


5.1 


8.1 


6 


196 


75.0 


0.52 


13.3 


22.6 


27 


200 


82.2 


0.56 


18.7 


31.0 


25 


202 


40.4 


0.48 


7.5 


13.0 


12 


205 


72.9 


0.54 


16.8 


27.9 


22 


207 


138.4 


0.52 


30.7 


51.9 


49 


Total 


639.0 


0.50 


108.5 


185.9 


165 


ZZ — > vvqq 
(Probabilistic selection) 




Integrated 


Selection 


Predicted 


Predicted 


Selected 


[GeV] 


luminosity [pb~^] 


efficiency 


background 


total MC 


data 


183 


52.4 


0.27 


1.8 


3.1 


4 


189 


152.8 


0.46 


24.3 


36.7 


32 


192 


24.9 


0.50 


5.6 


8.5 


9 


196 


75.0 


0.46 


12.7 


21.4 


20 


200 


82.2 


0.46 


16.4 


26.0 


25 


202 


40.4 


0.51 


10.9 


16.6 


16 


205 


72.9 


0.48 


18.0 


28.4 


24 


207 


138.4 


0.45 


29.6 


47.6 


48 


Total 


639.0 


0.45 


119.3 


188.3 


178 



Table 4: Integrated luminosities, selection efficiencies and number of observed and ex- 
pected events in the ZZ — > vvqq channel, after cuts on each discriminant variable max- 
imising the product of the efficiency and purity. In the IDA analysis, the measured 
ZZ production cross-sections were extracted at each centre-of-mass energy by fitting the 
distributions of the corresponding discriminant variable (see the text) . 
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variable, with the ZZ signal contribution as the only free parameter. The results obtained 
and the combinations performed between energies and with other channels to derive global 
values for the NC02 cross-section are described in Section 10. The optimal region of the 
distribution of the IDA variable to be used in the fit was determined by minimising 
the combined statistical and systematic uncertainty, evaluated using the full LEP-2 data 
sample. This resulted in removing events with values of the IDA variable less than 
2.2. The cross-section measured by combining the results of the probabilistic cross-check 
analysis over the different energies is also given in Section 10 for comparison, using the 
numbers of events and efficiencies presented in Table 4. 

The remaining background was composed of several processes. Near the threshold for 
ZZ production, the main backgrounds were from 55(7) events with badly reconstructed 
jets, and from 1^1^ processes. At higher energies, mis-reconstructed gg(7) events were 
easier to suppress, due to the Lorentz boost of the Z bosons in the signal. In this case, 
the main background was from single-resonant Wev^, processes. 

5.5 Systematic errors 

The main sources of systematic error in the selection of ZZ vvqq events were from 
uncertainties in the description of the energy flow reconstruction of the 55(7) background 
and, to a lesser extent, from uncertainties in the predicted background cross-sections, 
particularly those of the single-resonant Wev^, process. The propagation of these uncer- 
tainties to the final steps of the IDA analysis was studied. 

Large samples of events at the Z peak, collected in the same conditions as the high 
energy data, were used to evaluate the errors in the energy flow reconstruction. The selec- 
tion of ZZ — > vvqq events exploits the large missing energy characteristic of this channel, 
and hence is sensitive to the description of the low energy tail in the reconstruction of 
background processes such as qq{'^). The Z events were compared with the simulation 
to estimate corrections to the particle flow in several momentum bins, separately for the 
barrel and endcaps, and for charged and neutral particles. The corrections, consisting 
mainly of changes in multiplicities to account for observed efficiency losses and possi- 
ble duplication effects in the pattern recognition, were typically only a few per cent for 
charged particles, but could be larger in the case of neutral particles at low momentum 
or in the endcaps. It was found that the simulation overestimated particle reconstruction 
efficiencies in the highest momentum bins. This resulted in underestimating the lower 
tail of the total energy distribution. Application of the computed corrections signiflcantly 
improved the overall agreement between data and simulation at the pre-selection level for 
all energy flow observables such as the total charged and neutral energies, the visible mass 
or the acoplanarity. Partial improvement also resulted in the tails of the distributions 
most relevant to the selection of ZZ —>■ vvqq events, although comparisons were in this 
case limited by the statistics of the available samples. 

This correction procedure was then applied to all simulated high energy samples and 
the IDA analysis was repeated. The main effect was to increase the ^'^(7) background in 
the high purity regions relevant to the signal extraction by about one third on the aver- 
age, while other components were much less affected. The ZZ cross-sections were then 
fitted, at each ccntrc-of-mass energy, using these modified versions of the simulation. The 
differences in expected results were always of the same sign, corresponding to reductions 
in the ZZ cross-section obtained. After correcting for 1/2 of these discrepancies, ±1/2 
of them were used as conservative estimates of the errors resulting from uncertainties in 
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the energy reconstruction. The magnitude of the effects decreased with energy, because 
of the decreasing relative importance of the ^^(7) background. 

The impact of propagating these uncertainties in the energy reconstruction to the fitted 
ZZ — > qqvv cross-sections is shown in Tabic 5, using the full LEP-2 data sample. The 
combined effect from uncertainties in background cross-sections is also shown, assuming 
the following errors [21]: WW : ± 2%, Wev : ± 5%, 55(7) : ± 2%. Both types of error 
were assumed to be fully correlated between the energies analysed. 



[GeV] 


Systematic error from 
energy flow 


Systematic error from 

MC cross-sections 


183 


±14.5 % 


±4.0 % 


189 


±10.0 % 


±2.5 % 


192 


±7.5 % 


±2.5 % 


196 - 208 


±3.5 % 


±2.5 % 



Table 5: Expected effects on the ZZ — >■ uuqq cross-section fitted in the IDA analysis, from 
uncertainties in the energy flow reconstruction and predicted background cross-sections. 



6 Jets and a pair of isolated muons or electrons 

The iji^jji~qq and e^e'qq decay modes comprise 9.3% of the ZZ final states. The two 
final state leptons are typically well isolated from all other particles. This feature can 
be used to select such events with high efficiency and low background contamination in 
both the muon and electron channels. Events were selected initially without explicit cuts 
on the masses of the final state fermion pairs, in order to analyse simultaneously the 
production of ZZ and Z7* events and, in the case of the e'^e~qq channel, contributions 
from the (Z/7*)e"'"e~ process. Mass cuts were then apphed to isolate the ZZ contribution. 

6.1 Selection procedure 

A loose hadronic pre-selection was first applied, requiring that the events have at least 

7 charged particles and a charged energy above 0.30 ^/s. To reduce the radiative return 
to the Z boson, events were rejected if a photon with energy greater than 60 GeV was 
found. The selection procedures were then closely similar for both fj.^f^~qq and e^e~qq 
channels. In order to maximize the lepton identification efficiency, any charged particle 
with momentum exceeding 5 GeV/ c was considered as a possible lepton candidate, around 
which nearby photons, if present, could be clustered^. This was found necessary to 
improve its energy evaluation in the presence of final state radiation. In the case of the 
e^e'qq channel, photons with energy between 20 GeV and 60 GeV were also considered 
as possible electron candidates, to recover events in which the electron track was not 
reconstructed. 

Events with at least two lepton candidates of the same flavour, opposite charge (if 
none of the electron candidates originated from a photon) and invariant mass exceeding 
2 GeV/c^ were then selected. All particles except the lepton candidates were clustered 

"For this purpose any neutral particle identified as a photon with energy greater than 0.5 GeV was combined with 
the lepton candidate if the invariant mass of the resulting cluster did not exceed 0.4 GeV/c^. At most two photons were 
included in such a cluster, the photon giving the smallest mass increase being added first. 
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into jets using the JADE algorithm [31] (with y^i^ = 0.01) and a kinematic fit requiring 
four-momentum conservation was applied, correcting appropriately the errors on lepton 
energies in cases where photons had been added by the clustering procedure. 

At least one of the two lepton candidates was required to satisfy strong lepton iden- 
tification criteria, while weaker requirements were specified for the second. Muons were 
considered as strongly identified if selected by the standard DELPHI muon identification 
package [12], based mainly on finding associated hits in the muon chambers. For weak 
muon identification only a set of kinematic and calorimetric criteria was used. Electrons 
were considered as strongly identified when the energy deposited in the electromagnetic 
calorimeter exceeded 15 GeV or 60% of the total energy of particles composing the elec- 
tron candidate cluster and when the energy deposited in the hadron calorimeter was 
consistent with the expectations for electrons. For weak electron identification only re- 
quirements on the momentum of the charged particle in the cluster and on the energy 
deposited in the hadron calorimeter were used. Moreover electron candidates originat- 
ing from applying the clustering procedure around a photon were considered as weakly 
identified. 

Two discriminating variables were then used to select the events: p™"^ the lesser of 
the transverse momenta of the lepton candidates with respect to their nearest jet, and 
the per degree of freedom of the kinematic fit. 

6.2 Results 

The distribution of the mass of one fermion pair {Mu or Mgg) when the mass of the 
second pair is within 15 GeV/ c^ of the nominal Z mass is shown in Figure 7 (upper left and 
right-hand plots, respectively) for the selected events. The corresponding distribution of 
the sum of the masses of the two fermion-pairs is shown in the lower plot of Figure 7. The 
observed distributions arc in reasonable agreement with the predictions from simulation. 
To select on-shell ZZ production, cuts were placed simultaneously on the masses of the 
I'^l" pair, on the remaining hadron system and on their sum, taking into account in 
the performance optimization the different mass resolutions of these final states and the 
presence of the single-resonant process, e~^e~Z, in the case of the e^e~qq channel. The 
observed and predicted numbers of selected on-shcll ZZ events arc shown in Table 6. The 
residual background is composed of l^l^qq events generated outside the signal region (see 
Section 3), of other processes, principally W^W~, other ZZ decays and, in the case of 
e^e~qq, qq{'~f) production. These results were used to derive the combined values of the 
NC02 cross-section described in Section 10. 

6.3 Systematic errors 

Several possible sources of systematic errors were checked. Uncertainties in the ef- 
ficiency of the lepton identification were studied comparing semilcptonic WW events 
selected in data and simulation using the strong lepton identification criteria. Uncer- 
tainties in signal efficiencies from the description of the kinematic observables used were 
evaluated comparing the Pt and distributions in data and simulation for all l^l'qq 
events selected without mass cuts. Corresponding uncertainties in background levels were 
evaluated by comparing samples of events selected in data and in simulation, requiring 
both isolated tracks to be not identified as leptons, while maintaining all the other criteria. 
Finally, uncertainties in the background level in the e'^e'qq channel from fake electrons 
were studied with qq{'y) events selected in data and in simulation with purely kinematic 
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Table 6: Integrated luminosities, selection efficiencies and number of observed and ex- 
pected selected events in the ZZ ii'^ii~qq and ZZ e^e~qq channels. 
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Figure 7: Distribution of the mass of the lepton pair (hadron system) when that of the 
hadron system (lepton pair) is within 15 GeV/c^ of the nominal Z mass (upper left and 
right-hand plots, respectively) and of the sum of the masses of the lepton pair and hadron 
system (lower plot) for the selected events. The dots are the data taken in 1997-2000 
and the histograms represent the simulation. The background labelled l^l~qq {I = e, fi) 
corresponds to the predicted contributions from these final states not satisfying the signal 
definition described in Section 3. 



criteria. These effects and the statistical uncertainty of simulated data yielded combined 
systematic errors on the efficiency to select ^i~qq and e^e~qq events of ±3.0%, and on 
the relative uncertainty in the background level of ±15%^°. Finally, additional errors of 
±2% and ±3% were added in the procedure to extract and combine the cross-sections 
with other channels (see Section 10), to account for uncertainties from the signal defi- 
nition used (see Section 3), arising partly because of interference effects and because of 
biases resulting from the single-resonant (Z/7*)e^e^ process, respectively. 



7 Jets and a pair of isolated r-leptons 

The T^T~qq decay mode represents 4.7% of the ZZ final states. The search for T~^T~qq 
final states was based on the inclusive selection of isolated clusters of particles with low 

'^"These evaluations were limited in accuracy by the statistics of the available samples. 
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charged multiplicity and low invariant mass. Each pair of such clusters was considered as 
a possible r+r" candidate, with all other particles assumed to originate from a qq system. 
Assuming this topology, a set of discriminating variables reflecting the isolation of the 
clusters and the likelihood of the T^T~qq hypothesis was defined. All these variables were 
combined into a single quantity using a probabilistic method based on likelihood ratio 
products. Cutting on this combined variable enabled r^r'qq final states consistent with 
ZZ production to be selected with high efficiency and low background contamination. 

7.1 Definition of r-clusters 

Each charged particle with momentum exceeding 1 GeV/c was considered as a r-cluster 
candidate. Every other charged particle (with energy exceeding 1 GeV) and photon (with 
energy exceeding 0.5 GeV) was then successively added to each one of these candidates 
and the incremental change in the cluster mass AMd was computed. After sorting the 
particles in increasing order of AM^i, the particle which produced the minimal AM^i 
was kept if the charged and neutral particle cluster multiplicities were each less than 4, 
if AMci was less than 1.5 GeV/c^ and, for clusters with one (more than one) charged 
particle, if their mass, Md, was less than 2.5 (1.9) GeV/c^. 

After adding each particle, the determination of AMd was repeated for all remaining 
particles with the modified cluster. The procedure was continued until all combinations 
had been tried. If the initial charged particle was identified as an electron or a muon 
(using the standard packages described in [12]), only photons were added. To study 
the isolation of these r-cluster candidates, all the particles of the events were clustered 
into jets with the JADE algorithm [31] (using y^in = 0.01). The four-momenta of the 
tracks of the r-cluster were subtracted from the jets to which they belonged, and the 
nearest jet with a remaining energy exceeding 3 (1.5) GeV, for a cluster with one (more 
than one) charged particle, was determined. Two variables characterising the isolation of 
each cluster were defined: the angle 9d between the cluster and the nearest jet, and the 
transverse momentum of the most energetic particle in the cluster with respect to the 
nearest jet. 

7.2 Pre-selection 

Events were required to have at least 7 reconstructed charged particles and a charged 
energy greater than 0.30 ^/s. To suppress radiative returns to the Z, events were rejected 
if a photon with energy exceeding 55% of the beam energy was found. At least two 
r-clusters per event were required, satisfying the isolation criterion cos6d < 0.95 (0.92), 
for clusters with one (more than one) charged particle, having a sum of charges satisfying 
the condition [Qi + (^21 < 2, and with at least one of these clusters containing only one 
charged particle. 

All particles not used to define the two r candidates were then clustered into jets and 
a kinematic fit [25] with six constraints, consisting of four-momentum conservation and 
imposition of the condition M^+t- — Mgg — Mz, was applied to the event, where M^+^- 
and Mqq are the invariant masses of the r and jet systems, respectively. To exploit the 
large missing energy in the r decay, the missing energies of each r-cluster, (defined as 
the differences between each fitted cluster energy, E^^f^, and the corresponding measured 
value), were then required to satisfy: 

• min(L>^, D|) > -15 GeV, 
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• min(D]j, D^) > 10 GeV if the clusters contained identified leptons (electrons or 
muons) of the same type, 

• -D^ > 3 GeV (i = 1, 2) for r-clusters containing an identified lepton, 

• D]^ + > 10 GeV if one or both r-clusters contained an identified lepton. 

After all these pre-selection criteria, the total number of events observed in data and 
predicted in simulation were 505 and 460.7, respectively. 

7.3 Probabilistic selection 

A probabilistic method based on likelihood ratio products was used to combine several 
variables defined at the level of the r-clusters {Ej'^, D^e!^, cos 6*^/^, -P/'^), and the prob- 
abilities of three kinematic fits, performed using the constraints from four-momentum 
conservation, using in addition the hypothesis of ZZ production (as described in section 
7.2) or that of WW production (requiring the two di-jet masses reconstructed in the 
pairing with smallest to be equal to Mw)- The probability density functions were 
determined for each quantity by combining the predictions from all centre-of-mass ener- 
gies separately for r-clusters with leptons, with one charged particle and with more than 
one charged particle. The separation power of Y, the combined discriminant variable 
obtained, is illustrated in Figure 8, where the distribution of — log^^Q Y is shown for a 
sample enriched with signal events by requiring M^+^- and Mqq to be compatible with 
the Z boson mass. 



DELPHI 
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Figure 8: Distribution of the combined discriminating variable for the T^T~qq analysis, 
obtained after enriching with signal events by requiring Mr+,— and Mqq to be compatible 
with the Z boson mass. The dots are the data taken in 1997-2000 and the histograms 
represent the simulation. 
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7.4 Results 

To isolate the on-shell ZZ production process, a cut, chosen to maximise the product 
of the efficiency and purity, was apphed to the combined variable: — log^^Q^ > 0.5. The 
fermion pair mass distributions for the events selected by this cut (Mt-+t— and M^^) are 
shown in Figure 9, obtained after a kinematic ffi of the event with the constraints from 
four-momentum conservation only. The observed distributions are in reasonable agree- 
ment with the predictions from simulation. Additional mass cuts were used to improve 
the signal-to-background ratio for the final selection: 75 < M^+t— < 110 GeV/c^ and 
170 < M^+T— + Mqq < 200 GeV/c^. The observed and predicted numbers of selected 
events are shown in Table 7. The background is composed of qqqq, {e~^e~ , fi'^ fi~)qq and 
qq{'j) events. The qqqq final states were primarily from WW decays. The overall contri- 
bution from other ZZ decays was less than 30% of the total background. These results 
were used to derive the combined values of the NC02 cross-section described in Section 
10. 
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Figure 9: Distribution of Mt-+t~- and Mgg, obtained through kinematic fitting of the 
event with the constraints from four-momentum conservation only, for the selected events 
satisfying — logj^g ^ > 0-5 (see the text). The dots are the data taken in 1997-2000 and 
the histograms represent the simulation. 



7.5 Systematic errors 

The main sources of systematic error in the selection of ZZ r^r^qq events arose 
from uncertainties in the reconstruction of the r-clusters and associated isolation vari- 
ables. After the pre-selection, about 10% more events were selected in data as compared to 
Monte Carlo (see Section 7.2). At this level the selection is dominated by WW processes 
where one W decays into and the other into hadrons. Although the corresponding 
topology is not identical with that of the signal, it is similar since it contains a r-lepton 
to be identified in the presence of hadronic jets. A systematic uncertainty of ±10% on 
the selection efficiency was assigned to take this discrepancy into account and based on 
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Table 7: Integrated luminosities, selection efficiencies and number of observed and ex- 
pected selected events in the ZZ r^r^qq channel, after cuts on M^+^-, Mqq and on 
their sum. 

studies of uncertainties in the lepton identification used in the selection of the r-clusters. 
An uncertainty in the background level of ±30% was also specified because the fraction 
of background events originating from other ZZ decays (see Section 7.4) was neglected 
in the procedure to derive combined NC02 cross-sections. Both errors were assumed to 
be fully correlated between the centre-of-mass energies analysed. 

8 Final states with four charged leptons 

The decay modes with I — e^ji comprise 0.4% of the ZZ final states. The 

event topology is clean and the only significant background comes from non-resonant 
e^e~l'^l~ production. 

8.1 Selection procedure 

Events were first selected if they contained between 4 and 8 charged particles, accom- 
panied by at most 10 neutral particles, irrespective of particle identification. In order 
to take into account bremsstrahlung effects for candidate electrons, the momentum of a 
charged particle was replaced by the sum of energies of electromagnetic clusters in a nar- 
row cone around the track direction if it gave a larger value. The total invariant mass of 
the charged particles had to be greater than 50 GeV/c^, and the minimum invariant mass 
after discarding any one of the charged particles larger than 20 GeV/c^. All combinations 
of four charged particles with total charge zero were then examined, and a combination 
was selected if all the following criteria were satisfied: 

• all four tracks had impact parameters at the interaction point smaller than 3.0 and 
0.5 cm, respectively in the projections containing the beam axis and perpendicular 

to it, and polar angles between 10° and 170°, 

• at least three of the four charged particles had momenta greater than 5 GeV/c, and 
the least energetic particle had momentum greater than 2 GeV/c, 
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• a system of two oppositely charged particles was found with both their invariant and 
missing masses within 10 GeV/c^ of the Z boson mass and having the same flavour 
(when both were identified as either electrons or muons using the standard packages 
described in [12]), 

• the two particles complementary to this system were separated by at least 90° from 
each other, 

• the invariant mass of all pairs of oppositely charged particles in the event exceeded 
5 GeV/c^. 

8.2 Results 

The observed and predicted numbers of selected events are shown in Table 8. They 
are in overall agreement within the statistical errors. These results were used to derive 
the combined values of the NC02 cross-section described in Section 10. The two events 
selected at \fs = 188.6 and 199.5 GeV were in the ii'^ii~e^e~ and ii~ ii'^ ii~ final states, 
respectively. An event display for the latter is shown in Figure 10. 
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Table 8: Integrated luminosities, selection efficiencies and number of observed and ex- 
pected selected events in the ZZ — > l^l'l'^l" channel. 



8.3 Systematic errors 

The main sources of systematic error in the analysis of the ZZ — > I'^l'l'^l' decay 
modes with I = e, fi were from the statistics of the Monte Carlo samples produced to 
compute the residual background and the selection efficiency and from approximations 
used to translate the results in terms of NC02 cross-sections (see Section 3). These 
three uncertainties were propagated into the combination described in Section 10. They 
amounted to ±12%, ±5% and ±5% relative errors on the cross-sections obtained in this 
channel, respectively, and dominated over uncertainties in background cross-sections and 
reconstruction. The third error was assumed to be fully correlated between the energies 
analysed. 
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Figure 10: The ZZ /x+yU^yU+yU^ candidate selected at 199.5 GeV centre-of-mass energy. 
Three muons are observed in the barrel and one in the forward part of the detector (the 
muon track going to the left is out of the plane). A kinematic fit was performed assuming 
the production of two identical heavy objects, each decaying into [i^ [T . Out of two 
possible combinations to form the fi~^fi~ pairs the highest probability was found for the 
top-bottom and left-right pairs. The corresponding fitted invariant mass, 90.8 GeV/c^, 
is fully compatible with ZZ production. 



9 Two isolated electrons or muons with missing en- 
ergy 

The vvl^l^ decay modes with I = e, fi comprise 2.7% of the ZZ final states. The 
events are characterized by two acoUinear charged leptons of the same flavour, with both 
invariant and recoil masses close to the Z mass, and by large missing energy. Although 
it has different production kinematics, the WW process also contributes to these final 
states with a large cross-section. A significant fraction of the corresponding events have 
exactly the same features and constitute a substantial background. 

9.1 Selection procedure 

To ensure good reconstruction, tracks were required to have impact parameters at the 
interaction point smaller than 3.0 and 0.5 cm, respectively in the projections containing 
the beam axis and perpendicular to it, and polar angles between 20° and 160°. Events 
with two particles identified as e'^e" or /i (using the looser identification criteria 
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described in [12] and calorimetric requirements), were selected if their total energy was 
less than 60% of y^, if the angle between them was in the range [^^j„,170°], where 
^min — 0.95arccos(l — 8m|/s), if the polar angle of the event missing momentum was 
between 25° and 155°, and if the reconstructed invariant masses satisfied: 

min{| Mz - m{l+l-) \ , \ Mz - rrirecaiiil^l') |} < 4 GeV/c^, 
max{| Mz - m{l+l-) \ , \ Mz - mrecoiill'^l') |} < 8 GeV/c^ 

where mrecoii{l^l ) is the invariant mass recoiling against the pair. Final state 
radiation for candidate electrons was taken into account as discussed in Section 8.1 for 
the four-lepton analysis. 

9.2 Results 

The observed and predicted numbers of selected events are shown in Table 9. They 
are in overall agreement within the statistical errors. These results were used to derive 
the combined values of the NC02 cross-section described in Section 10. 
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Table 9: Integrated luminosities, selection efficiencies and number of observed and ex- 
pected selected events in the ZZ — > i/P/i+//~, vDe'^e~ channels. 



9.3 Systematic errors 

The main sources of systematic error in the analysis of the ZZ — > uul'^l' decay modes 
with I = e, fi were from the statistics of the Monte Carlo samples produced to compute 
the remaining background and the selection efficiency and from uncertainties in the back- 
ground level. These three uncertainties were propagated into the combination described 
in Section 10. They amounted to ±13%, ±5% and ±5% relative errors on the cross- 
sections obtained in this channel, respectively, and dominated over experimental errors 
in the selection efficiency and over uncertainties from approximations used to translate 
the results in terms of NC02 cross-sections. The third error was assumed to be fully 
correlated between the energies analysed. 
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10 Combined NC02 cross-sections 

The measurements described in the previous sections for the different ZZ channels 
show reasonable agreement between the data and simulation. In this section the compar- 
ison is presented in terms of extracted NC02 cross-sections combined over all channels at 
each centre-of-mass energy, as well as over all energies by normalising to the SM predic- 
tions, both separately in each channel, and for all channels together. Such combinations 
enable more comprehensive and statistically meaningful checks of the SM prediction. 
Cross-sections for individual sub-channels are also given at each energy in the case of the 
most important ZZ final states. 

10.1 Likelihood method 

To obtain the cross-section information in the different cases, probability functions at 
each energy and for each sub-channel, defined with respect to variations of the value of 
the NC02 cross-section, were combined into global likehhoods. For the channels qqqq and 
vvqq (analysed in the IDA stream), the shapes of the probability functions were derived 
from the fits to the distributions of the combined variables defined in each analysis. 
For the remaining channels the Poissonian probabilities were constructed, based on the 
numbers of events selected in the data and predicted in the simulation. The central 
value for each measurement was defined as the median of the corresponding likelihood. 
The statistical errors were obtained from the two intervals on either side of the median 
each containing 34.13% of the total probability. In cases when less than 31.74% of the 
probability was below the maximum of the likelihood, an upper limit was given, defined 
as the value which was exceeded by 5% of the integrated probability distribution. 

10.2 Propagation of systematic errors 

The propagation of systematic uncertainties affecting each different final state to its 
probability function was studied by introducing appropriately modified assumptions for 
backgrounds and efficiencies (as described in the corresponding sections). The following 
procedure was applied to propagate all systematic errors in the combination of any selec- 
tion of channels or energy points. Uncertainties were divided into four groups according 
to their correlation among energies and channels. The only significant channel-correlated 
uncertainties were on the integrated luminosities collected at each centre-of-mass energy. 
They amounted to ±0.3% (theory error) and ±0.6% (measurement error) and were re- 
spectively treated as fully correlated and uncorrelated among the energies. 

A set of Gaussian random numbers was then drawn and used to assign values for each 
kind of error taking into account its nature (common random number among channels 
or energies in cases of fully correlated errors, and independent ones otherwise). The 
probability functions were then modified accordingly and their combination was repeated 
a large number of times. The standard deviation of the Gaussian-like distribution of 
central values obtained was taken to represent the total systematic error affecting the 
combined measurement for the selection of channels or energy points considered. 

10.3 Results and discussion 

The measured values of the NC02 cross-section at each centre-of-mass energy are hsted 
in Table 10 and shown in Figure 11, together with the SM predictions from the YFSZZ 
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Figure 11: Combined NC02 cross-sections measured from data collected in 1997-2000. 
The errors shown are sums in quadrature of the statistical and systematic errors. The solid 
curve is the SM prediction computed using the YFSZZ generator [32]. The corresponding 
uncertainty is about ±2% [21]. Values obtained with the EXCALIBUR generator were 
consistent with these predictions within 1%. 



generator [32]. These measurements, normalised to the predicted SM values, were also 
combined into a single value: 

Rnco2 = 0.91 ± 0.08 (stat) ± 0.02 (syst). 

The systematic errors at individual energies and in the combination had the same mag- 
nitude because the main uncertainties were fully correlated. The measurements at ^/s = 
183 and 189 GeV were compatible with those previously obtained [7]. The uncertainty in 
the SM prediction, of about ±2% [21], is not shown (as for the other results in this sec- 
tion). A good overall agreement with the SM prediction can be observed. The expected 
weights of each measurement in the combination are given in Table 10. 

The combination of the results from the different energies was also performed for each 
channel separately. The results are listed in Table 11, with their expected weight in 
the global combination, and with the corresponding ZZ branching fractions predicted 
in the SM. The results are also shown graphically in Figure 12. The channels were 
ordered according to their expected weight. This order was not always correlated to 
the decay branching ratios because of the different background conditions. As can be 
seen, the three main channels ZZ qqqq, vi'qq and l~^l~qq (with I = e, fi) had similar 
weights, and accounted for over 90% of the total sensitivity. A good overall agreement 
with the SM predictions can be observed. The combination was also performed for the 
probabilistic uuqq cross-check analysis, and the result obtained: Rnco2 = 0.80^Q;}g (where 
the errors are statistical), was consistent with that obtained from the main IDA analysis. 
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Table 10: Measured and predicted NC02 cross-sections at the different energy points. In 
the column with measured values, the first errors are statistical and the second system- 
atic. The SM prediction was computed using the YFSZZ generator [32]. The uncertainty 
on these predictions is about ±2% [21]. Values obtained with the EXCALIBUR generator 
were consistent with these predictions within 1%. The last column contains the expected 
weight of each energy point in the global combination performed over all energies de- 
scribed in the text. 

This channel had the largest systematic uncertainty of all (relative to statistical errors). 
However, in this channel as in all others, systematic uncertainties were much smaller than 
the statistical errors. In the four-lepton channel, only an upper limit (quoted at 95% CL) 
is given. The of the seven measurements had a probability of 70%. 

Finally, for completeness, the values of the measured cross-sections normalised to the 
SM predictions are also presented at each energy in Table 12 for the three dominant 
channels. 
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Table 11: Ratios of measured to predicted cross-sections for individual channels combined 
over all energy points (for the two last channels / = e,/i). The predicted cross-sections 
used are listed in Tabic 10. The upper limit given for the channel corresponds to 

a 95% confidence level. In the column with measured values, the first errors are statistical 
and the second systematic. The two last columns contain the weight of the channel in 
the global combination over all channels described in the text, and the corresponding ZZ 
SM branching fraction. 
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Table 12: Measured cross-sections normahsed to the SM predictions for individual chan- 
nels at individual centre-of-mass energies. The predicted cross-sections used are listed in 
Table 10. The quoted errors are statistical. The upper limits shown correspond to a 95% 
confidence level. 
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Measured/Predicted Cross-Section 

Figure 12: Ratios of measured to predicted cross-sections for individual channels com- 
bined over all energy points (for the two last channels / = e,/i). The predicted 
cross-sections used are listed in Table 10. The errors shown are sums in quadrature 
of the statistical and systematic errors. The vertical band displays the total error on 
the combination of the seven channels. The upper limit given for the channel 
corresponds to a 95% confidence level. 
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